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ABSTRACT 

New black-hole mass estimates are presented for a sample of 72 AGN covering three 
decades in optical luminosity. Using a sub-sample of Seyfert galaxies, which have black- 
hole mass estimates from both reverberation mapping and stellar velocity dispersions, 
we investigate the geometry of the AGN broad-linc region (BLR) . It is demonstrated 
that a model in which the orbits of the line-emitting material have a flattened geom- 
etry is favoured over randomly orientated orbits. Using this model we investigate the 
Mth— Lbuige relation for a combined 90-object sample, consisting of the AGN plus a 
sample of 18 nearby inactive elliptical galaxies with dynamical black- hole mass mea- 
surements. It is found that, for all reasonable mass-to-light ratios, the M^h — Lbuige 
relation is equivalent to a linear scaling between bulge and black-hole mass. The best- 
fitting normalization of the Mth — M^uige relation is found to be M^h — 0.0012Mfc„;ge, 
in agreement with recent black-hole mass studies based on stellar velocity dispersions. 
Furthermore, the scatter around the M^h — Lbuige relation for the full sample is found 
to be significantly smaller than has been previously reported (AlogMb/j = 0.39 dex). 
Finally, using the nearby inactive elliptical galaxy sample alone, it is shown that the 
scatter in the M^h — Lbuige relation is only 0.33 dex, comparable to that of the M^h — c 
relation. These results indicate that reliable black-hole mass estimates can be obtained 
for high redshift galaxies. 

Key words: galaxies: active - galaxies: nuclei ~ galaxies: bulges - quasars: general 
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1 INTRODUCTION 

The technique of using the broad H/3 emission line to esti- 
mate the central black-hole masses of active galactic nuclei 
(AGN) has recently been employed widely in the literature 
(eg. Lacy et al. 2001, Laor 2001a, Wandel 1999). Given the 
ease with which the nuclear spectra of AGN can be obtained, 
a proper calibration of line-width black-hole mass estimates 
has the potential to allow the study of the evolution and de- 
mographics of active supermassive black holes over a wide 
range in redshift. 

In a previous paper (McLure & Dunlop 2001, hereafter 
MDOl) we reported H/3 black-hole mass estimates for a 45- 
object sample consisting of both Seyfert galaxies and pow- 
erful {Mv < —23) quasars. By combining the black-hole 
mass estimates with host-galaxy bulge luminosities derived 
from full two-dimensional disc/bulge decompositions, MDOl 
found that quasars and Seyfert galaxies follow the same 
A^bh — Ltuige relation, a result recently confirmed by Laor 
(2001a). Furthermore, by adopting a simple inclination cor- 
rection factor for the H/3 line widths, MDOl found that the 



mean Mbh/Mbuige ratio in AGN is a factor of 2 — 4 lower 
than had been determined by Magorrian et al. (1998) for 
nearby galaxies, in reasonable agreement with the results 
from recent stellar velocity dispersion studies (eg. Merritt & 
Ferrarese 2000a). 

Despite the recent attention which has been focussed 
on determining the black-hole masses of both active and 
inactive galaxies, several important problems remain to be 
resolved, all of which are potentially soluble. 

Firstly, at present the usefulness of the Mth ~ Lbuige 
relation as a black-hole mass estimator for both active and 
inactive galaxies is severely limited due to its large scat- 
ter (~ 0.5 dex). In our previous study (MDOl) we demon- 
strated that much of this scatter could be due to the dif- 
ficulty of accurately determining the bulge luminosities of 
late-type galaxies, even at 2; < 0.1. Although the correla- 
tion between black-hole mass and stellar velocity dispersion 
for nearby inactive galaxies displays a much smaller scatter 
(~ 0.3 dex, Merritt & Ferrarese 2000b), it is still clear that 
a Mbh — Lbuige correlation with reduced scatter would be 
highly desirable, given the extreme difficulty in obtaining 
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stellar velocity dispersions for high redshift galaxies. In this 
paper we investigate the possibility that the intrinsic scat- 
ter in the Mhh — Lbuige relation is substantially lower than 
previously reported by studying a sample of nearby inac- 
tive galaxies for which Lbuige has been determined to high 
accuracy. 

A further complication which arises when studying the 
Albh — Lbuige relation for powerful active galaxies, is the un- 
certainty in how best to calibrate the virial black-hole mass 
estimates produced by using the broad H/3 line-width to 
trace the central gravitational potential. Indeed, uncertainty 
about the exact geometry of the broad-line region (BLR) in 
AGN has the potential to produce large systematic errors in 
line-width based black-hole mass estimates (see Krolik 2000 
for a discussion) . In this paper we directly address this issue 
by adopting a flattened disc-like BLR geometry, which we 
demonstrate to be fully consistent with presently available 
data. 

Thirdly, it is currently unclear whether the Mbh — Mbuige 
relation is linear over a large baseline in bulge mass. Al- 
though our previous study (MDOl) found no evidence for 
non-linearity in the Mbh — Mbuige relation, the recent study 
of Laor (2001a) came to the opposite conclusion, finding 
that Mbh oc MfJ^fgg. As well as being of intrinsic interest, 
the question of the linearity of the Mbh — Mbuige relation can 
also be used to constrain current models of combined black 
hole/bulge formation (eg. Kauffmann & Haehnelt 2000). 

Fourthly, the results of our HST host-galaxy study 
(McLure et al. 1999, Dunlop et al. 2001) together with 
our previous H/3 study (MDOl), and the II/3 study of Laor 
(2001b), point to an apparent division between the black- 
hole masses of optically selected radio-loud and radio-quiet 
quasars at ~ 10^ M© . However, recent results from Lacy 
et al. (2001) and Dunlop et al. (2001) both indicate that 
this apparent threshold is not of fundamental physical im- 
portance. In this paper we use a large sample of 72 AGN, 
together with a sample of 20 nearby inactive galaxies with 
dynamic black-hole mass estimates, to systematically re- 
address these questions. 

The structure of this paper is as follows. In Section ^ 
we describe the construction of the AGN and inactive galaxy 
samples. In Section |^ we briefly review the technique of esti- 
mating black-hole masses from emission-line widths, before 
proceeding to describe our adopted disc BLR model, and 
then to demonstrate that it is fully consistent with recent 
velocity dispersion measurements for a subset of our Seyfert 
galaxy sample. In Section ^ the normalization and linearity 
of the resulting Mbh — Lbuige relation for the AGN and inac- 
tive galaxy samples is investigated and compared with litera- 
ture results. In this section we also compare the Mbh — Lbuige 
and Mbh — relations for the inactive galaxy sample, and 
demonstrate that the scatter in both is virtually identical. In 
Section ^ we investigate the correlation between black-hole 
mass and the ratio of bolometric luminosity to the Edding- 
ton limit for the AGN sample. In Section H we explore the 
implications of the new black-hole mass estimates for the ra- 
dio loudness dichotomy, before presenting our collusions in 
Section ^. To provide consistency with our previous study, 
all cosmological calculations in this paper assume Ho = 50 
kms-^Mpc-\ n = 1.0, A = 0. 



Sample 


N 


bu I g e 


H/3 FWHM 


QSO 


30 


MDOl 


MDOl 


QSO 


11 


H97 


FOl 


QSO 


8 


MPDOl 


BG92 


QSO 


4 


POl 


BG92 


Syl 


15 


MDOl 


WPM99 


Syl 


4 


VOO 


LOl 


Inactive 


18 


F97 





Table 1. Details of the data drawn from the literature. Column 2 
lists the number of objects in each sub-sample. Column 3 lists the 
principle references for the bulge luminosity data. Column 4 lists 
the references for the H/3 line-width data. The reference codes 
are as follows: MDOl (McLure & Dunlop 2001), H97 (Hooper et 
al. 1997), FOl (Forster et al. 2001), MPDOl (McLure, Percival 
& Dunlop 2001), BG92 (Boroson & Green 1992), POl (Percival 
et al. 2001), WPM99 (Wandcl, Peterson & Malkan 1999), VOO 
(Virani et al. 2000), LOl (Laor 2001a), F97 (Faber et al. 1997). 



2 THE SAMPLE 

The full sample analysed in this paper consists of 90 ob- 
jects comprising three sub-samples of 53 quasars, 19 Seyfert 
1 galaxies and 20 inactive nearby galaxies. This sample is 
a combination of the 45 objects from our previous study 
(MDOl), together with 47 additional objects drawn from 
various literature sources. The principle objective behind 
the construction of the sample was to allow the study of the 
bulge: black-hole mass relationship over the widest possible 
dynamic range in both nuclear luminosity and central black- 
hole mass. Due to the fact that the sample is drawn from 
several sources, it was necessary to transform all the bulge 
luminosity data into the same filter. As in MDOl, the cho- 
sen filter was the standard _R-Cousins, which was preferred 
because measured _R— band bulge luminosities were avail- 
able for 65/72 of the AGN sample. The principle references 
for the bulge luminosity and line-width data are provided 
in Table |l|, with speciflc details relating to each sub-sample 
provided below. 

2.1 The quasar sample 

The 53 objects in the quasar sample cover the redshift range 
0.1 < 2 < 0.5, and have absolute magnitudes ranging from 
Mr ~ -22.5 to Mr ~ -28.0. Consequently, the lowest 
luminosity quasars in the sample cover the overlap region 
around the Seyfert /quasar divide, while the highest lumi- 
nosity objects constitute some of the most powerful quasars 
available at 2; < 0.5. Thirty of the objects are drawn from 
the quasar sample analysed in MDOl, which have accurate 
7?— band bulge luminosities from the host-galaxy study of 
Dunlop et al. (2001), and have line- width and optical con- 
tinuum measurements from either our own recent observa- 
tions (MDOl), or the study of Boroson & Green (1992). 
The remaining 23 quasars in the new sample are drawn 
from three additional sources. Eleven further objects are 
taken from the HST host-galaxy study of LBQS quasars by 
Hooper et al (1997). The bulge luminosities are converted 
from the Rj magnitudes published by Hooper et al. assum- 
ing Rc — R.J = 0.1 (Fukugita et al. 1995), with the line- 
width and continuum measurements taken from Forster et 
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al. (2001). Data for four objects is taken from the host- 
galaxy study of high-luminosity radio-quiet quasars by Per- 
cival et al. (2001). These quasars have A"— band bulge lumi- 
nosities based on sub-arcsecond tip/tilt imaging obtained at 
UKIRT, and optical spectra taken from Boroson & Green 
(1992). The final eight objects are taken from a new study 
(McLure, Percival & Dunlop 2001, in prep) of the host- 
galaxies of a radio-loud companion sample to that studied 
by Percival et al. The bulge luminosities in this study are 
taken from a combination of new A"— band UKIRT observa- 
tions and archive HST data. The optical spectra for these 
objects are also taken from Boroson & Green (1992). The 
bulge luminosities for those objects with Ji'— band data were 
converted to the _R— band assuming a rest-frame colour of 
7? - if = 2.5 (Dunlop et al. 2001). 

2.2 The Seyfert galaxy sample 

The Seyfert galaxy sample features the 15 objects analysed 
in MDOl, with the addition of four objects from the Seyfert 
galaxy study of Virani et al. (2000) for which line- width data 
was available. In addition to the four new objects, the Rc 
disc/bulge decompositions provided by Virani et al. have 
also been adopted for three objects from the MDOl sample 
(NGC 4051, NGC 4151, and NGC 3227). The bulge lumi- 
nosities for these objects from Virani et al. replace those 
originally used by MDOl, which were taken from Baggett, 
Baggett & Anderson (1998), and are preferred here because 
the data were obtained using the Rc filter. The bulge lu- 
minosities for the other 12 objects remain the same as in 
MDOl. 



2.3 The inactive galaxy sample 

The normal galaxy sample consists of eighteen objects drawn 
from the list of 37 nearby inactive galaxies with dynamical 
black-hole mass measurements published by Kormendy & 
Gebhardt (2001). Due to the difficulty in determining accu- 
rate bulge luminosities in late-type galaxies (see MDOl for 
a discussion) it was decided to exclude those galaxies in the 
Kormendy & Gebhardt list which were not of E-type mor- 
phology (including lenticulars) . This decision was taken to 
investigate to what extent the scatter in the Mbh — Lbuige 
relation could be reduced if the uncertainties in the bulge 
luminosities were minimized (see Section ^ . The Kormendy 
& Gebhardt list consists of 20 E-type galaxies, which was re- 
duced to the final 18-object sample by the exclusion of NGC 
4486B and NGC 5845. These two objects were excluded be- 
cause the errors on their black-hole mass measurements are 
an order of magnitude larger than for the rest of the sample 
(we note that NGC 4486B was also excluded as an outlier by 
Merritt & Ferrarese (2000b) in their study of black-hole de- 
mographics). The _B— band bulge luminosities for 9 of the 18 
objects were taken from Faber et al. (1997), with 1/— band 
data for a further 9 being taken from Merritt & Ferrarese 
(2000b). To convert the _B— band and 1/— band magnitudes 
to the i?— band, standard bulge colours of B — R = 1.57 
and V — R = 0.61 were assumed (Fukugita et al. 1995). All 
eighteen objects in the sample have published velocity dis- 
persion measurements, allowing a direct comparison of the 
Mbh — Lbuige and Adbh ~ <y relations (see Section ^). 



3 THE VIRIAL BLACK-HOLE MASS 
ESTIMATE 

The technique of using H/3 line-widths to trace the gravita- 
tional potential of the central black holes which power AGN 
has been used extensively in recent years. Detailed discus- 
sions are provided by Wandel, Peterson & Malkan (1999) 
and Krolik (2001). The basic premise is that the dominant 
mechanism responsible for the width of the broad emission 
lines is the gravitational potential of the central black hole, 
and that the line-widths reflect the keplerian velocities of 
the line-emitting material (see Peterson & Wandel (2000) 
for supporting evidence). If this assumption is valid then 
the so-called virial mass estimate for the central black hole 
is given by: 

Mbh = RblrV'G-^ (1) 

where Rblr is the radius of the BLR and V is the veloc- 
ity of the line-emitting material. As in MDOl, the method 
adopted in this paper for estimating the BLR radius is the 
correlation between Rblr and monochromatic luminosity 
at 5100A found by Kaspi et al. (2000), from a combination 
of the reverberation mapping results for 17 PG quasars and 
17 Seyfert galaxies (15 of which are in the sample studied 
here). Incorporating this into Eqn|l| yields the following for- 
mula for the black-hole mass: 

Mbh = 3.64 X f X °" X {FWHMf (2) 

where Mbh is the black-hole mass in solar units, FWHM is 
the full- width half maximum of the H/3 line in kms~^, and / 
is a inclination factor which links the observed H/3 FWHM 
to the intrinsic velocity of the line-emitting material. It can 
immediately be seen from Eqn|^that the inclination factor / 
can have a potentially large effect on the derived black-hole 
masses. Although evidence exists that radio- loud AGN dis- 
play inclination dependent FWHMs consistent with a disc- 
like BLR (see below), in the radio-quiet population the ge- 
ometry of the BLR is arguably undetermined. Consequently, 
the standard assumption made in the literature is that the 
velocities of the line-emitting material are randomly orien- 
tated, which leads to / = (eg. Wandel 1999). However, 
as in our previous study (MDOl), in this paper we adopt 
a specific model for the geometry of the broad-line region, 
which, as discussed below, is more consistent with the avail- 
able data. 

3.1 The disc model of the broad-line region 

Following MDOl we make the assumption in this paper that 
the broad-line emitting material has a flattened disc-like ge- 
ometry, and consequently that the observed H/3 FWHM de- 
pends on the orientation of the disc rotation axis relative 
to the line of sight. Substantial evidence exists in the lit- 
erature that the FWHM of broad emission lines are incli- 
nation dependent in radio-loud AGN (eg. Wills & Browne 
1986, Brotherton 1996, Corbin 1997, Vestergaard, Wilkes 
& Barthel 2000). Such strong observational evidence does 
not exist for radio-quiet AGN, largely because the absence 
of radio jets removes one of the main methods with which 
to constrain the AGN orientation axis. However, the exten- 
sion of a disc BLR model to the radio-quiet AGN in our 
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V„t/1000kms ' FWHM/lOOOkms ' 

Figure 1. The left-hand panel shows the location (cross) of the best-fitting values of Vorh and o-^rb^ together with the 1,2&3(T joint 
confidence levels defined by Ax^ = 2.3, 6.2 & 11.8. Also shown (filled circle) is the location of the best-fitting parameters obtained by 
minimizing the Kolmogorov-Smirnov statistic between the model and data cumulative FWHM distributions. The right-hand panel shows 
the inclination correction factor (/) derived from the best-fitting parameters (solid line), together with the itlo" models (dashed lines). 
Also shown is the effective inclination factor resulting from the assumption of random broad-line velocities (horizontal dotted line) and 
the mean observed FWHM of our sample (vertical dotted line). 



sample can be justified on several grounds. Firstly, under 
the unified scheme (eg. Urry & Padovani 1995) it is ex- 
pected that (to first order) the central engines, and presum- 
ably BLR, of radio-loud and radio-quiet AGN are essentially 
identical. Secondly, as will be discussed further below, pre- 
vious studies have demonstrated that it is possible to fit the 
broad emission-line spectra of both radio-loud and radio- 
quiet AGN using a BLR model with a fiattened geometry. 
Thirdly, as will be seen in Section 3.3, for those radio-quiet 



objects in our sample where it is possible to test the geome- 
try of the BLR, it appears that the data are most consistent 
with a fiattened geometry. 

In order to determine the parameters of the chosen disc 
model we simply require that the model can reproduce the 
distribution of the observed H/3 FWHMs. In MDOl, the cu- 
mulative FWHM distribution of a smaller 45-object sample 
was fitted using a two-population model, where both sub- 
populations had the same characteristic orbital velocity, but 
were confined to lie within difi'erent ranges of inclination to 
the line of sight. Although the model adopted in MDOl pro- 
vided an excellent fit to the available data, it was not fully 
satisfactory. The main reason for this is that because there 
was no way to individually allocate objects to one of the 
two populations, it was therefore only possible to calculate 
an average inclination correction, based on the sample as a 
whole. Consequently, it was therefore inevitable that objects 
with FWHM at the extremes of the observed distribution 
would not be properly corrected for inclination, systemati- 
cally increasing the scatter in the resulting black-hole mass 
estimates. 

Consequently, in the analysis of this new larger 72- 
object AGN sample it was decided to investigate the pos- 
sibility of providing an acceptable description of the sam- 
ple FWHM distribution using a single population model, 
uniquely determined by only three free parameters. This 
model makes the assumption that the intrinsic keplerian or- 
bital velocities of the sample have a gaussian distribution, 
with mean Vorb and variance a^^b- The third free parameter 



between the line of sight and the disc rotation axis. In the 
model fit it is presumed that the angles between the line 
of sight and the disc rotation axis are randomly distributed 
between zero (pole-on) and 6max- In order to relate the in- 
trinsic orbital velocities to the observed FWHMs we follow 
the results of the accretion disc modelling of Rokaki & Bois- 
son (1999). Rokaki & Boisson used an accretion disc model 
to fit the UV continuum and H/3 line profiles of a sample of 
19 Seyfert 1 galaxies (9 of which are common to our sample) 
with the fitting process returning estimates of both central 
black-hole mass and the disc inclination angle. From their 
results for the 9 objects in common to both samples it fol- 
lows that: 



FWHM ~ 2 X sinie) [GMbh/Rsat 



(3) 



of the model is 



which is the maximum allowed angle 



where Rsat is the radius at which the radial disc line-fiux 
Fi{r) saturates. It was shown by Rokaki & Boisson that 
the derived values of Rsat are in good agreement with 
the broad-line radii Rblr derived from reverberation map- 
ping observations. In this study we therefore make the as- 
sumption that Rblr = Rsat, which immediately leads to 
FWHM = 2Vorb X sin{9). Consequently, in this disc model 
the inclination correction factor is simply / — 2sit,{e) • 



3.2 Modelling the FWHM distribution 

The parameters of the best-fitting model are determined by 
minimizing the chi-squared statistic between the model and 
data FWHM distributions. The best-fit parameter values 
(x?nin = 2.75, 5 d.o.f.) determined from this process are 
as follows: Vorb = 4375 kms~^, aorb ~ 1400 kms~^ and 
Omax = 47°. The location of the minimum chi-squared so- 
lution, together with the relevant confidence level contours, 
is shown in the left-hand panel of Fig ^ As a consistency 
check we have also determined the best-fit parameters by 
minimizing the Kolmogorov-Smirnov (KS) statistic between 
the model and data cumulative FWHM distributions. The 
location of this minimum is also indicated in Fig ^, and can 
be seen to be consistent with the minimum chi-squared so- 
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cr/kms ^ Reference 



FWHM/lOOOkms 



-1 



Figure 2. The cumulative FWHM distributions of the AGN sam- 
ple (thick line) and the best-fitting disc BLR model (thin line). 
The two can be seen to be indistinguishable (KS, p = 0.99) 



lution. The quality of the model fit can be seen from Fig ^ 
which shows the cumulative FWHM distribution of the full 
72-object sample and the best-fitting model. The applica- 
tion of the KS test confirms that the two distributions are 
essentially identical {p = 0.99). 

The second panel of Fig ^ shows the inclination fac- 
tor / predicted by the best-fitting disc model, together with 
the predictions of the ±lcr models. Also shown in the plot 
(horizontal dashed line) is the effective inclination factor re- 
sulting from the assumption of randomly orientated BLR 
orbits; / = The prediction of random orientated BLR 
orbits has been included to illustrate that the disc BLR 
model adopted here only makes significantly different pre- 
dictions to the random case for those objects with observed 
FWHM< 2800 kms"\ For objects with FWHM> 2800 
kms~^ (83% of the 72-object AGN sample) the inclination 
correction factor derived from our disc model differs by less 
than a factor of 2 from that resulting from the assumption 
of random velocities. 



3.3 Testing the geometry of the broad-Hne region 

Although the inclination corrections introduced by the disc 
BLR model only have a significant effect for a minority of 
objects in the AGN sample, it will be demonstrated in the 
next section that they significantly reduce the scatter in the 
- Li^ige relation. 

Due to this fact, it is important that the inclination de- 
pendence of the observed FWHM predicted by the disc BLR 
model can be shown to be at least consistent with the avail- 
able data. As was discussed in Section 3.1, although there 
is substantial evidence in the literature for an inclination 
effect in the observed FWHM of radio-loud AGN, such evi- 
dence is not yet available for radio-quiet AGN. As a result, 
the justification for extending the disc model to radio-quiet 
objects relies primarily on the assumption that the central 
engine of all powerful AGN are essentially the same. 

However, included in our 19-object Seyfert galaxy sam- 
ple are 10 objects for which reliable stellar velocity disper- 
sion measurements are available in the literature (see Table 
2). Because stellar velocity dispersions are orientation inde- 



3C 120 


162 


NW95 


MRK 79 


120 


FOl 


MRK 590 


169 


NW95 


MRK 817 


140 


FOl 


NGC 3227 


128 


NW95 


NGC 3516 


144 


A97 


NGC 4051 


88 


NW95 


NGC 4151 


85 


FOl 


NGC 5548 


180 


FOl 


NGC 6814* 


115 


NW95 



Table 2. The ten Seyfert galaxies in our sample which have pub- 
lished stellar velocity dispersions. All of the objects, with the ex- 
ception of NGC 6814, also have reverberation mapping measure- 
ments of Rblr- For NGC 6814 Rblr has been estimated from 
its continuum luminosity. Column three lists the literature sources 
for the velocity dispersions, with the following codes: NW95 (Nel- 
son & Whittle 1995), FOl (Ferrarese et al. 2001) and A97 (Arribas 
et al. 1997). 



pendent, it follows that a direct comparison of the black- hole 
masses predicted by the Mbh — o" relation with those of the 
orientation dependent H/3 method will allow a test of the 
BLR geometry. 

The procedure is as follows. For each of the ten Seyfert 
galaxies we first estimate their black-hole mass using the 
best-fitting M^h — o" relation to our nearby inactive galaxy 
sample (see Section W) which has the form: 



logMbh = 4.09(±0.58) log(a) - 1.26(±1.39) 



(4) 



a relations 
and 



and can be seen to be consistent with the M^h 
determined by both Merritt & Ferrarese {M^h oc a 
Gebhardt et al. (2000b)(Mi,h oc ct^'^'^). We then require for 
each object that this stellar velocity dispersion black-hole 
mass estimate be equal to the H/3 virial black-hole mass 
estimate given by Eqn |l|, using the reverberation mapping 
estimate of Rblr and assuming V = /xFWHM. This pro- 
vides a model-independent estimate of the inclination fac- 
tor / for each of the 10-objects. In Fig ^ we plot the es- 
timated values of / against observed FWHM for the ten 
Seyfert galaxies, along with the predicted curve from the 
disc BLR model, and the effective inclination factor result- 
ing from the assumption of randomly orientated BLR veloc- 
ities. From Fig ^ there can be seen to be good agreement 
between the derived inclination correction factors and those 
predicted by our adopted disc BLR model, particularly in 
the right-hand panel where the data have been binned. Also 
shown in the right-hand panel of Fig ^ is the inclination 
correction curve produced by fitting to the FWHM distribu- 
tion of Seyfert galaxies alone. Due to its lower characteristic 
black-hole mass, it can be seen that this curve represents an 
improved match to the trend displayed by the binned data. 
We note here that the possible impact of fitting the var- 
ious AGN sub-samples with different characteristic orbital 
velocities has been investigated. However, because the dif- 
ferences introduced into the black-hole mass estimates are 
small (~ 0.1 dex), this more complicated procedure has neg- 
ligible effect on the results presented in the rest of the paper. 
Consequently it was decided to adopt the best-fitting model 
to the full AGN sample throughout. 
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FWHM/lOOOkms ' FWHM/lOOOkms ' 

Figure 3. The left-hand panel shows the derived inclination factor (/) versus FWHM for the 10 Seyfert galaxies with measured stellar 
velocity dispersions (see text for discussion). Also shown is the prediction of the best-fitting disc BLR model to the full AGN sample 
(solid curve) and the expected value of / in the case of purely random BLR velocities (dotted horizontal line). The error bars have been 
estimated by assuming a scatter of 0.3 dex in the Mi,h — a relation and ±15% errors in FWHM measurements. The right-hand panel 
shows the same information using binning to clarify the trend displayed by the data. The dashed curve shows the inclination prediction 
of the best-fitting disc model to the complete 12/im Seyfert 1 sample of Rush, Malkan & Spinoglio (1993), which has a lower mean 
black- hole mass than our full AGN sample (see text for discussion). 



Although the small sample size makes it unwise to draw 
any firm conclusions, we simply note that the assumption 
that the observed FWHM of the AGN sample are inclination 
dependent is at least consistent with the available data, and 
apparently more so than the assumption of purely random 
BLR velocities. 



4 THE RELATIONSHIP BETWEEN BULGE 
LUMINOSITY AND BLACK-HOLE MASS. 

In Fig ^ absolute J?— band bulge magnitude is plotted 
against black-hole mass for the 72 objects in the AGN sam- 
ple. Also shown is absolute 7?— band bulge magnitude plot- 
ted against dynamically-estimated black-hole mass for the 
20 objects in our nearby inactive galaxy sample. Two as- 
pects of Fig ^ are worthy of immediate comment. Firstly, 
as was shown by MDOl and Laor (1998 & 2001a), it can 
be seen that bulge luminosity and black-hole mass are ex- 
tremely well correlated, with the Spearman rank correla- 
tion test returning Vs = —0.77 (7.3(t). Secondly, it is clear 
that the AGN and nearby inactive galaxy samples follow 
the same M^h — Ltuige relation over > 3 decades in black- 
hole mass, and > 2.5 decades in bulge luminosity. This sec- 
ond fact strongly supports the conclusions of Dunlop et al. 
(2001) and Wisotzki et al. (2001), that the host-galaxies of 
powerful quasars are normal massive ellipticals drawn from 
the bright end of the elliptical galaxy luminosity function. 
Thirdly, there can be seen to be no systematic offset between 
the Seyfert 1 and quasar samples, reinforcing the findings of 
MDOl and Laor (2001a) that, contrary to the results of Wan- 
del (1999), the bulges of Seyfert galaxies and QSOs form a 
continuous sequence which ranges from (bulge) ~ —18 
to M_R(bulge)~ —24.5. If we adopt an integrated value of 

= -22.2 (Lin et al. 1998), then this implies that the 
Mbh — Lbutge relation holds from Lbutge — O.OIL*, all the 
way up to objects which constitute some of the most mas- 
sive ellipticals ever formed; Lf,uige — lOL*. 

In order to find the best-fitting Mth — Ltuige relation 



we use the same iterative method used in MDOl, which 
allows for measurement errors in both coordinates (Press et 
al. 1992). The best-fit to the full 90-object sample has the 
following form: 

log(A//i,h/ Mq) = -0.50(±0.02)Mii - 2.96(±0.48) (5) 

and is shown as the solid line in Fig ^. The scatter around 
this best-fitting relation is l^M^h = 0.39 dex. It is worth 
noting that this level of scatter means that the Mbh — Ltuige 
relation is now worthy of increased interest, given that it 
amounts to an uncertainty factor of < 2.5. Given that the 
20 objects in the nearby inactive galaxy sample have actual 
dynamical black-hole mass estimates, it is obviously of in- 
terest to quantitatively test how consistent the Mth — Ltuige 
relation for these objects is with the fit to the full, AGN 
dominated, sample. The best-fit to the inactive galaxy sub- 
sample alone, has the following form: 

log(Mi,h/ Mo) = -0.50(±0.05)Mr - 2.91(±1.23) (6) 

which can be seen to be perfectly consistent with the best-fit 
to the full sample in terms of both slope and normalization. 
Indeed, as shown in Table ^ the best-fitting relations for 
the full sample, quasar sample, Seyfert galaxy sample and 
the nearby inactive galaxy sample are all internally con- 
sistent, and display comparable levels of scatter. This is a 
remarkable result given that it implies that the combined 
bulge/black hole formation process was essentially the same 
throughout the full sample, which as well as featuring both 
active and inactive galaxies, includes galaxies of both late 
and early-type morphology. 

The quality of the fit to the inactive galaxy sample is 
illustrated by the left-hand panel of Fig |^, which shows the 
Mth — Ltuige relation for the inactive galaxy sample alone. 
Of particular interest is the scatter around this best-fit re- 
lation, given that it has been widely reported in the lit- 
erature (eg. Merritt & Ferrarese 2000b, Kormendy & Geb- 
hardt 2001) that the scatter around the Mth — Ltuige relation 
is significantly greater than that around the Alth — cr rela- 
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Sample 


N 


a 


b 


9 


AA4h 


All 


90 


-2.96 ±0.48 


-0.50 ±0.02 


1.43 


0.39 


Active 


72 


-2.25 ±0.72 


-0.46 ± 0.03 


0.94 


0.42 


QSOs 


53 


-4.26 ±3.15 


-0.55 ±0.13 


0.87 


0.42 


Seyfert 


19 


-1.96 ± 1.33 


-0.45 ±0.06 


1.21 


0.43 


Inactive 


20 


-2.91 ± 1.04 


-0.50 ±0.05 


1.05 


0.33 



Table 3. Details of the best-fitting relations of the form 
log Afj,h = a + hM[i{bulge) to the full sample and to the various 
sub-samples. Column two gives the number of objects included in 
each fit. Columns three and four give the best-fitting values of the 
parameters a ii b, along with their associated uncertainties. The 
relatively large uncertainties in the fit to the QSO sub-sample re- 
flects the small dynamic range in terms of both -LbuZge 

and Mj,h. 

Column five gives the reduced foi' each fit and column 6 gives 
the scatter around each of the best-fitting relations in terms of 
AlogM^^. 

tion. However, in contrast, we find that the scatter around 
the Mi,h — Ltuige relation for our sample of nearby inactive 
galaxies, which excludes non E-type morphologies, is only 
0.33 dex, in excellent agreement with the scatter around the 
Mbh - o relation (Merritt & Ferrarese 2001b). 

To test this result further, in the right-hand panel of 
Fig ^, we investigate Mbh — o" relation for our nearby inac- 
tive galaxy sample. The scatter around the best-fit relation 
{Mhh oc cr^"^) is 0.30 dex, leading us to the conclusion that 
the intrinsic scatter around the M^h — Ltuige relation for el- 
liptical galaxies is comparable to that in the M^h — o" relation. 
This may largely be due to the fact that when construct- 
ing the nearby inactive galaxy sample, objects with late- 
type morphologies were deliberately excluded. This would 
support the conclusion of MDOl that successful disc/bulge 
decomposition of late-type galaxies is a difficult task, even 
with high resolution data, and that bulge luminosities from 
ground-based images with poor seeing and, in particular, 
morphology-based estimates of bulge/total luminosity frac- 
tion (eg. Simien & de Vaucouleurs 1986) can often have sub- 
stantial systematic errors associated with them. However, 
we also note that the late-type galaxies excluded from our 
nearby galaxy sample are not significant outliers with re- 
spect to the Mth — o" relation. Consequently, although our 
results indicate that bulge luminosity can provide accurate 
black-hole mass estimates, comparable with the Mbh — o" re- 
lation, for elliptical galaxies, it is clear that the Mi,h—<^ rela- 
tion will provided more accurate black-hole mass estimates 
for samples which are not restricted to solely ellipticals. 



4.1 The linearity of the bulge: black-hole mass 
relation 

As pointed out in Section one of the aims of this paper 
was to investigate the linearity of the Mth ~ M^uige relation. 
In our previous study (MDOl) of a sample of 45 AGN we 
found that Mhh oc Afj^jg^!^'''^^, and therefore concluded that 
there weis no evidence that the Mhh — Mhuige relation was 
non-linear. In contrast, evidence for a non-linear relation 
was recently found by Laor (2001a). In his band study 
of the black hole to bulge mass relation in a 40-object sample 

© 0000 RAS, MNRAS 000, 000-000 



CD 





1 ' 1 

o 






1 












A o 


o 


C\? 


■a ns jf 


°o ° 




cv 

1 


■ 




- 

- 


o 










■ . ■ 






1 




A 


Normal galaxies 






■ 


Seyfert Is 


CO 




o 


QSOs 


1 


— ■ y " 

/I , 1 




1 




6 8 




10 



Figure 4. Absolute _R— band bulge magnitude versus black-hole 
mass for the full 90-object sample. The black-hole masses for the 
72 AGN are derived from their H/3 line-widths under the disc- 
like BLR model. The black-hole masses of the inactive galaxies 
(triangles) are dynamical estimates as compiled by Kormendy & 
Gebhardt (2001). Also shown is the formal best-fit (solid line) 
and the best-fitting linear relation (dotted line). 

(15 PG quasars, 16 inactive galaxies and 9 Seyfert galaxies) 
Laor found a best-fitting relation of the form: 

Mhh=MtTgf-'' (7) 

which is clearly apparently inconsistent with linearity. How- 
ever, in order to determine the Mi,h — M^uige relation it 
is obviously necessary to convert the measured bulge lu- 
minosities into masses, via an adopted mass-to-light ratio. 
The form of this mass-to-light ratio affects the derived slope 
of the Mhh — Mbuige relation in the following way. If the 
mass-to- light ratio is parameterized as M/L cx L", then the 
resulting slope (7) of the Mth — Mhuige relation is given by 
7 — , where /3 is the slope of the Mhh — Lhuige relation 

(Eqn 5)? 

As in MDOl we choose to adopt the derived ii— band 
mass-to-light ratio for the Coma cluster from j0rgensen, 
Franx & Kjaergaard (1996), which has a = 0.31. With 
this mass-to-light ratio the best-fitting Mbh — Lbuige rela- 
tion (Eqn ^ transforms to a Mbh — Mbuige relation of the 
following form: 

MbhOcMttf°-°' (8) 

It can immediately be seen that that there is no indication 
from our results that the relation is inconsistent with a linear 
scaling between black hole and bulge mass. 

In order to calculate the bulge mass of the objects in 
his sample, Laor (2001a) adopted a 1/— band mass-to-light 
ratio of Mbuige oc L\'^ig^ (Magorrian et al. 1998), which is 
significantly different from our chosen mass-to- light ratio. 
Indeed, if Laor (2001a) had adopted the mass-to-light ratio 
used here, then his best-fit to the Mbh — Mbuige relation for 
his full sample would be Mbh oc M^^^^^^ ' , which is not 
formally inconsistent with linearity. However, irrespective of 
this, our new best-fit to the slope of the Mbh — Lbuige re- 
lation (/3 = —0.50 ± 0.02) of our new sample, which has a 
larger dynamic range in Lbuige than both the samples stud- 
ied in MDOl and Laor (2001a), means that any disagreement 
about mass-to-light ratios cannot now alter the conclusion 
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Figure 5. Left-hand panel shows absolute R— band bulge magnitude versus dynamical black-hole mass estimate for our 18-object inactive 
galaxy sample. The solid line is the best-fitting relation which has the form M^jt oc j^j^^g^J^""^ ■ The dotted line is the best-fitting linear 
relation to the full sample, and is equivalent to Af^/i = 0.0012Mj,ujge- The right-hand panel is the same with bulge luminosity replaced 
by stellar velocity dispersion. The solid line is the best-fit (M;,;, cx: cr* '^^) , the dashed line is the Merritt &; Ferrarese (2000b) relation 
{Mf,ii oc (T*'^^), and the dot-dashed line is the Gebhardt et al. (2000b) relation (M^jj cx: a^-"^^). In both figures the location of the Milky 
Way and M31 are indicated for the interest of the reader, although neither was included in the nearby galaxy sample for the purposes 
of the analysis. 



that the Mth — Mi,uige relation is consistent with being lin- 
ear. To demonstrate this we conclude by noting that even 
using the Mj,„;ge oc L\'^^g^ mass-to-light ratio adopted by 
Laor (2001a), our best-fitting Af^h — Ltuige relation is equiv- 
alent to: 

Mbh OC Mj„,g^ (9) 
again, completely consistent with a linear scaling. 

4.2 The normalization of the Mbh — Mbuige 
relation 

Having established that the M^h — Mi,aige relation is con- 
sistent with being linear, we now assume perfect linearity 
in order to establish the normalization of the Mth — Mtuige 
relation. With the mass-to-light ratio adopted here, a lin- 
ear scaling corresponds to enforcing a slope of —0.524 in the 
Mbh vs. Mr relation. Under this restriction the best-fitting 
relation (dotted line in Figs ^ & 5a) has a normalization of 
Mbh ~ 0.0012Mbuige, and can clearly be seen to be an ex- 
cellent representation of the data. It is noteworthy that the 
normalization of Mbh ~ 0.0012Mi,„ij5 is identical to that de- 
termined by Merritt & Ferrarese (2000a) from their velocity 
dispersion study of the 32 inactive galaxies in the Magorrian 
et al. sample. 

The closeness of the agreement between the Mbh /Mbuige 
ratios determined here with those determined by Mer- 
ritt & Ferrarese is highlighted by Fig ^ which shows 
a histogram of the Mbh /Mbuige distribution for our 72- 
object AGN sample. The AGN Mbh /Mbuige distribution 
has {\og{Mbh/ Mbuige)) = -2.87 ± 0.06 with a standard 
deviation oi a = 0.47. This is in remarkably good agree- 
ment with the Merritt & Ferrarese results, which were 
{\0Z{Mbh/ Mbuige)) = -2.90 and a = 0.45. 

The close agreement between the distribution of 
Mbh /Mbuige found here for the 72-object AGN sample and 
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Figure 6. Histogram of the ratio of black-hole mass to bulge 
mass for the 72-object AGN sample. Over-plotted for compari- 
son is a gaussian with {log{Mbh/ Mbuige)) = —2.90 and standard 
deviation 0.45 (see text for discussion) 

that found by Merritt & Ferrarese using velocity dispersions 
of nearby inactive galaxies, can also be taken as further sup- 
porting evidence that the adoption of a disc model for the 
BLR geometry is valid. Finally, we note that the normaliza- 
tion of Mbh = 0.0012Mf,uige agrees very well with the predic- 
tions of recent models of coupled bulge/black hole formation 
at high redshift (Archibald et al. 2001). 



5 THE BOLOMETRIC LUMINOSITY - 
BLACK-HOLE MASS CORRELATION 

In this section we use the optical continuum information 
available for the AGN sample to investigate how bolomet- 
ric luminosity as a fraction of the Eddington limit varies 
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Figure 7. Black-hole mass plotted against bolometric luminosity 
as a fraction of the Eddington limit for the 72-object AGN sample. 
Radio-loud quasars are shown as filled circles, radio-quiet quasars 
are shown as open circles and the Seyfert galaxies are shown as 
open squares. The best-fitting relation is shown as the solid line, 
with the location of the Eddington limit shown by the vertical 
dotted line. Also shown is a representative error bar. 



with black-hole mass. In order to estimate the bolomet- 
ric luminosity of each AGN we follow Wandel, Peterson & 
Malkan (1998) and adopt Lboi — lOALsioo, where L5100 is 
the monochromatic luminosity at 5100 A. Although this cal- 
ibration was also found to be appropriate for the composite 
quasar spectrum of Laor et al. (1997), it is obviously only 
a rough estimate (Wandel, Peterson & Malkan 1998) and 
consequently we allow for a factor of 4 uncertainty when 
deriving the best-fitting relation. 

In Fig ^ we plot black-hole mass against Lboi/ LEdd for 
the full 72-object AGN sample. It can be seen that the two 
quantities are well correlated, with a Spearman rank coeffi- 
cient of Ts = 0.50, 4.2(T. The best-fitting relation (solid line) 
has the following form; 

logMth = 2.79(±0.37) \og{Ltoi/ L,,da) - 10.29(±0.28) (10) 

and shows clearly that the data are inconsistent with all 
AGN having a constant Eddington ratio, a conclusion which 
was also arrived at by Kaspi et al. (2000). As in Fig ^ it 
is apparent that the Seyfert galaxies and quasars form a 
continuous sequence, and that in terms of Eddington ra- 
tio, there is no evidence for any systematic offset between 
the radio-loud and radio-quiet quasars. However, it is note- 
worthy that there is no evidence for a correlation within 
the QSO sub-sample alone. Indeed, given that the observed 
correlation is entirely dependent upon the heterogeneous 
Seyfert galaxy sample, it should perhaps be viewed with 
some caution. For example, it is inevitable that objects 
which should occupy the top-left corner of Fig ^ will be 
missing from our sample, since their large bulge luminosity 
and weak nuclear emission will have prevented them from 
being classified as powerful AGN. However, it is not immedi- 
ately obvious why there should be no members of the Seyfert 
galaxy sample which are radiating at much closer to their 
Eddington limit. Although it is possible that the bottom- 
right of Fig could be populated by so-called narrow-line 
Seyfert galaxies, which are thought to have intrinsically 
small black-hole masses and to radiate at a large fraction of 



their Eddington limit (eg. Mathur 2000), our present study 
does not support this. Four of the objects in the Seyfert 
galaxy sample have FWHM< 2200 kms~^ and yet do not 
occupy this region of the diagram. The implications of the 
Mbh — Lboi /LEdd relation, together with the Mbh-Mbuige re- 
lation determined in the previous section, will be explored in 
a forthcoming paper (McLure, Percival & Dunlop, in prep). 



6 BLACK-HOLE MASS AND THE 
RADIO-LOUDNESS DICHOTOMY 

In this final section we investigate the relationship between 
black-hole mass and radio luminosity for the objects in the 
AGN sample, to explore the role played by black-hole mass 
in the quasar radio-loudness dichotomy. To ensure that any 
connection between black-hole mass and radio power can be 
investigated properly it is essential that the radio-quiet and 
radio-loud sub-samples should be indistinguishable in terms 
of both redshift and optical nuclear luminosity. From the 
total of 53 quasars in the AGN sample it has been possible 
to construct two sub-samples of radio-loud and radio-quiet 
quasars, numbering 23 and 26 objects respectively, which 
are well matched in term of redshift and optical luminosity 
(ALsioo). The matching of the two sub-samples is confirmed 
by the application of the two-dimensional KS test (Peacock 
1983), which returns a probability of p = 0.25 that the two 
distributions are indistinguishable in the AL5100 — z plane. 
Given that the quasars as a whole are a heterogeneous sam- 
ple, this AL5100 — z matching is the best that can be done 
to ensure that any differences in the sub-sample black-hole 
mass distributions should be linked to the large differences 
in radio power. 

The mean black-hole mass of the radio-quiet quasar 
sample is (log(Mf,,,/ M©)) = 8.61 ±0.07, with a median 
value of log(Mf,(,/ Mq) = 8.56. The black-hole masses of 
the radio-loud quasar sub-sample are somewhat larger on 
average, with a mean of (log(A'/t.ji/ M0)) = 8.76 ± 0.07, 
and a median of log(M(,h/ M©) — 8.83. The suggestion from 
these results is that there does appear to be a tendency for 
the radio-loud quasars to have larger black-hole masses than 
their radio-quiet counterparts (the median figures differ by 
nearly a factor of 2), although the overlap between the two 
sub-samples is sufficient to ensure that the two black-hole 
mass distributions are not statistically distinguishable (KS, 
p — 0.21). Given that the radio-loud and radio-quiet sub- 
samples were deliberately chosen to be matched in terms 
of their nuclear optical luminosity, it is clear from Eqn ^ 
that the difference in average black-hole mass between the 
two sub-samples must be due to differences in observed H/3 
FWHM. This does appear to be the case, and can be seen 
in Fig ^, which shows a plot of optical luminosity against 
observed FWHM for the two sub-samples. As is highlighted 
in Fig ^ there is a separation around FWHM=5500 kms~^, 
with 11/23 RLQs having FWHM > 5500 kms"\ while only 
4/26 RQQs have FWHM> 5500 kms'^ This difference is 
confirmed as being significant, with the KS test returning a 
probability of p = 0.05 that the two FWHM distributions 
are drawn from the same parent distribution. The difference 
in the FWHM distributions between radio-loud and radio- 
quiet quasars has been previously reported at higher signif- 
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icance {p = 0.001) from the larger-scale spectrophotometric 
studies of Corbin (1997) and Boroson & Green (1992). 

The suggestion that the black holes of radio-loud 
quasars are on average more massive than their radio-quiet 
counterparts is consistent with the results of our recent host- 
galaxy study (Dunlop et al. 2001, McLure et al. 1999), our 
previous H/3 line-width study (MDOl) and the recent line- 
width study of the LBQS by Laor (2001b). In this study, 
Laor calculated virial black-hole mass estimates for the 87 
LBQS objects with z < 0.5, finding an apparent separation 
between the radio-loud and radio-quiet quEisars at a black- 
hole mass of ~ 10® M©. However, it is important to remem- 
ber that the radio-loud objects in both of these studies are 
predominantly optically selected. The possible weakness of 
using optically selected quasar samples to study the link be- 
tween black-hole mass and radio power has recently been 
highlighted by Lacy et al. (2001). In their line-width study 
of a virtually complete sub-sample of the FBQS, Lacy et 
al. find that the apparent gap in the radio: black-hole mass 
plane is filled in by radio intermediate quasars, and that the 
quasax population as a whole follows a relation of the form: 

Although the black-hole masses and 5GHz radio lumi- 
nosities of the radio-loud quasars analysed here are consis- 
tent with the results of Lacy et al., the large scatter around 
the Lacy et al. relation (1.1 dex) makes this somewhat in- 
evitable. In fact, we argue elsewhere (Dunlop et al. (2001), 
McLure & Dunlop, in prep) that the data are more consis- 
tent with the existence of an upper and lower envelope to 
the radio power that can be produced by a black hole of 
a given mass. Consequently, we argue that in order to pro- 
duce a truly radio-loud AGN {Psghz ^ 10^^ WKz-^sr-^), 
the central black-hole mass nmst be > 10*^"'' Mq and lie on 
the upper Mi,h — P^ghz envelope (see Dunlop et al. (2001) 
for a full discussion) . Finally, we note that this picture of the 
radio-loudness dichotomy is in good agreement with the pre- 
dictions of the merger model proposed by Wilson & Colbert 
(1995), in which the radio-loudness of an AGN is dependent 
on both the mass and angular momentum of the central 
black hole. 



7 CONCLUSIONS 

New virial black-hole mass estimates have been presented 
for a sample of 72 AGN using the FWHM of the broad H/3 
emission line, corrected for inclination under the assumption 
of a disc-like broad-line region. The form of the Mi,h — Mtuige 
relation is investigated using reliable bulge luminosity mea- 
surements for both the AGN sample and a carefully selected 
sample of neaxby inactive galaxies with dynamical black-hole 
meiss measurements. Using the inactive nearby galaxy sam- 
ple a rigorous comparison was made of the scatter around 
the Mih — Lbuige and M^h — a relations. The main conclu- 
sions of this study can be summarized as follows: 

• For a 10-object subset of the Seyfert galaxy sample it is 
found that the inclination predictions of a disc model BLR 
are more consistent with the available data than the assump- 
tion of purely random broad-line velocities. 

• A strong correlation (7.3(t) is found between bulge lu- 
minosity and black-hole mass estimated via the disc BLR 
model. 
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Figure 8. The distribution of the radio-loud and radio-quiet 
quasar sub-samples in tlie AL5100— FWHM plane. It should be 
noted that there is a significant diflference in the FWHM distri- 
butions of the two samples (KS, p = 0.05), with 11/23 radio-loud 
quasars having FWHM> 5500 kms~^ while only 4/26 radio-quiet 
quasars have FWHM> 5500 kms"!. 



• The best-fitting Mth — Ltuige relation to the combined 
sample of 72 AGN and 18 nearby inactive elliptical galaxies 
is found to be consistent with a linear scaling between black 
hole and bulge mass 

lower scatter than previously reported (Alog Mbh = 0.39 
dex). 

• The best-fitting normalization of the Mbh — Mbuige rela- 
tion is found to be Mbh = 0.0012Mbuige, in excellent agree- 
ment with recent stellar velocity dispersion studies. 

• In contrast to previous reports it is found that the scat- 
ter around the Mbh — Lbuige and Mbh — c" relations for the 
nearby inactive elliptical galaxy sample are comparable, at 
only ~ 0.3 dex. 

• Using samples matched in terms of optical luminosity, 
the median black-hole mass of radio-loud quasars is found 
to be larger by 0.27 dex than that of their radio-quiet coun- 
terparts. This difference is found to be due to a small but 
significant (p = 0.05) difference in their respective FWHM 
distributions. 
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